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Abstract 

The approach used for the determination of S-wave amphtudes containing the apphcation 
of the nonrelativistic approximation in the case of P-waves leads to unphysical divergen- 
cies. We show how to avoid the latter in calculations of contributions to the cross section 
near threshold in agreement with field theory. This enables us to give quantitatively re- 
liable predictions for the forward-backward asymmetry and for the axial contribution to 
the total cross section for the top-antitop system. Also the cross section for the production 
of stop-antistop near threshold is determined. 
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1 Introduction 



The observation of the top quark has been reported recently by two colaborations at 
FNAL 0, g (CDF and DO). They quote masses of 176 ± 8 ± lOGeV and 199^^? ± 22GeV 
respectively . Electroweak data from LEP and the predicted Standard Model cross 
section indicate a slightly lower top mass around 170GeV. 

Thus it seems to be established by now that the top quark is very heavy, giving it an 
extraordinary status in the particle spectrum discovered so far. A main characteristic of 
the top quark compared to other heavy quarks is its large decay width due to t ^ W^b 
of the order of IGeV. In the present paper we will investigate properties of the top which 
are relevant in a future e~^e~ collider (or perhaps fi^fi~). 

Such a machine would be the ideal place to investigate in detail its properties. Es- 
pecially the threshold region could give us independent determinations of the mass and 
width of the top quark as well as an improved determination of the strong coupling con- 
stant. Whereas for lighter quarks the qq bound state shows individual energy levels (as 
in charmonium and bottonium) the large electroweak decay width of the top smears out 
individual levels. On the other hand it makes QCD perturbation theory more reliable 
01 because the t decays fast enough on the times scale of strong interaction to prevent 
hadronisation and thus essential contributions of inherent nonperturbative character. 

This work is organized as follows. In sect. 2 we fix the notation for cross sections and 
summarize the Green function approach. Section 3 is devoted to the calculation of the 
forward-backward asymmetry originating in the interference of S- and P-waves. P-wave 
Green functions are needed exclusively for the calculation of the total cross section for 
the stop-antistop production as well as for the axial contribution to the cross section for 
ti. Both are calculated in chapter 4. Finally in chapter 5 we draw our conclusions. 
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2 Calculation of cross sections near threshold 



The nowadays quite standard Green function approach [§, ^ rehes on the observation 
that particles produced near threshold move relatively slow and thus one approximates 
the exact four point function by a nonrelativistic Green function. This Green function 
fulfills the equation 

{-— + V{r) -E- iT}G{f, f') = -5{f- f') (1) 
m 

where V{r) is a perturbative or QCD-motivated potential. The width F has been included 
to treat an unstable particle. In general the total cross section for ti production is given 
by ( a color factor 3 already included, X = A,V ) 

a = ay + aA (2) 
72 

ax = -pj7r[cx + rfxp]ImG^^(P)a^^, (3) 

whereas for scalar constituents, like tt (stop-antistop) from an eventual supersymmetric 
generalization only one contribution occurs 



72 

= p2 t^'' + dsp]l^Gs{P)(T^p (4) 



with 

:= (aW)2 + (6W)2 = + b 



(5) 

dx := 2a(^)6(^) 4 = 2aA 



and 



P = 2m + E. (6) 



The Standard Model (SM) values of a^^\b^^^ as well as the values of the minimally 
supersymmetric extended Standard Model (MSSM) as,bs are given by 

VtVeP'^ , (V) VtaeP'^ 



_ atVeP^ _ ataeP^ . . 
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and 



where we used the abbreviations 

z = P"^ - Ml + iMzTz, s = sin^vF, c = cosew (10) 
Vf = Ti-2Qjsm^ew, af = Ti (11) 

Here tI is the eigenvalue of the diagonal SU(2) generator for the fermion / (e.g. t or 
e)with charge Qf. For the supersymmetric charges we use |p 

= ~qu (12) 

= {cosset -2qtsm^ew)/{2sc). (13) 

The leading contributions to (0) arises from the vector coupling since it produces ti 
pairs with angular momentum zero (S-waves). In the fermionic case P- waves are produced 
by the axial current and represent higher order correction. They contribute to 0{a1) to 
the total cross section and lead to a forward backward asjmametry by interfering 
with the vector contribution. On the contrary, for scalar particles P-waves represent the 
dominant contribution near threshold. The four point functions G^^{P),Gs{P) are in 
the nonrelativistic approximation replaced by 



ImG^^(P) ^ -^lmG{0, 0) (14) 

toG-(P)^i_M^|_G(J,,/)|,„^„,„_„ (15) 

1 . d d ^ 



ImGsiP) ^ ;;^^"^^^^(^'^)l|-Ho,|,i^o (16) 

While ImG(0, 0) in eq. (|l^ is well behaved, the r.h.s in eqs. ([15|) and ([16|) are linearly 
divergent. 

Befor tackling this problem, let us consider the calculation of the nonrelativistic Green 
function. Due to its symmetry and the requirement of regularity at r — oo and at the 
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origin, the general solution of eq. (Q) can be written as 

oo I 

Gir,r') = Y.9iir,r') Yr^in')YU^) (17) 

/=0 m.=-l 

I n 9<{r<)g>{r>) \ r:r>r' \ r : r < r' 

gi{r,r) = , r> = <^ r< = i (18) 

where g<^{r^) and g^{ry) are regular solutions of the homogeneous equation 



- m{E + zr - V{r)) - ^-^^}9{r) = (19) 



at r = or r ^ oo, respectively. The actual behavior of the solution g<{r^) (and also of 
the irregular one) for r ^ depends on the potential and on /. In this work we use a one 
loop renormalization group improved potetial, treating four flavors as massles (nj = 4). 
The bottom quark contribtion is included in a pure perturbative manner. 

TT/ X cP' 5 5 1 , , 5,, , , 

= r..-2n.w. . , ^ + —\E<-rm,e-^ ' 7: + + e^)] (20) 



r[l-^2^^^(7 + ln/ir)] ^vrr^ ^ 6 2 

The color factor 4/3 is absorbed in the definition of a: 



b 



4 o2 

« ■■= = ^ (21) 



In principle all other parts of the general potential derived in |TT| on purely field theoretic 
grounds can be included as well. 

Due to the presence of the width iV the behavior of (7>(r) for r — > cxd is given by 

lim g>ir) oc e"""^ (22) 

with 

a_ := ^^-E + + r2. (23) 

This exponential damping behavior is the origin of the infrared cut off mentioned already 
in IQ. Consider for simplicity the case E' = 0. Then we can write 



where Eb = and tb are the Bohr energy and the Bohr radius of the system, 

respectively . This formula clarifies that if the width F is approximately equal to Eb the 
form of the potential is only "tested" up to the order of magnitude of the Bohr radius. 
The situation is improved for E < —T since here we have 



V m\E\ 

This relation also holds also for small F and thus the Bohr radius again becomes the 
relevant scale. 

The matching condition for and is provided by the 6 distribution in eq. ([l|): 

-m = g<{r)g'^{r) - g'^{r)gy{r) (25) 

Numerically it is possible to obtain the regular solution at the origin directly by imposing 
suitable boundary conditions. In the following we will need both S- and P-wave Green 
functions. For the singular solution of the former we use the following method p, |^. 
Suppose you have two solutions: (7<(r) as above and U2{r) a solution determined by an 
arbitrary boundary condition. Then the solution g^ is given by 

^>(r) = cK(r)+5(7<(r)] (26) 

with 

5 = - hm (27) 



The constant c can be determined from eq. (pSj). For I = and a Coulomb- like potential 
at the origin we arrive with the boundary conditions 

g^{0) = (28) 
9<iO) = 1 

and arbitrary ones for U2 at the result 

Attw r 

w = U2g'^ - g<U2. (30) 
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According to eq. (|^) and (|T3D the total cross section is proportional to ImG(0,0). With 
this knowledge it is possible to calculate the total cross section near threshold for any 
given potential V{r). We have written a numerical routine for MATHEMATICA"^^^ which 
has been checked to give the correct answer for the known Green function of a purely 
Coulombic potential. 

Before presenting some numerical studies, a short remark on the coupling constant to 
be used seems in order. It is customary to take the strong coupling constant as in the MS 
scheme. Especially experimental determinations are always given in terms of ag^^{Mz)- 
But in bound state calculations it is natural to use an as defined differently. Clearly it 
should be possible to relate the two schemes. To avoid complications from heavy fermion 
masses we treat 4 flavors as massless. The difference in the bottom quark contribution 
can be expected to be smaller than the experimental uncertainties in as- Therefore we 
ignore it in the following consideration. The quark-antiquark potential in the MS scheme 
can be extracted e.g. from [l^ to 0{as): 



Ana ,^ , 31 lOrif, 

V = -^ l-a(37raoln^ + — -— ^) 

Comparison with the QCD potential |]11| for nj light flavors gives 



(31) 



«"^ = «M5 (1-^(31-^)) (32) 



where BS denotes our bound state scheme used in |[TT|]. Numerically this means that our 
a is slightly smaller than usual which is advantageous for our perturbative calculation. 
For rrit = 180, a^^{Mz) = .117 ± 0.05 get for a renormalization at /i = I/tb^h) 

/i = lS.h±lGeV (33) 
a*^(/i) = 0.22 ±0.01 (34) 
a(/i) = 0.19 ±0.01 (35) 

Remember that the above values for a, a^^ differ by a factor 4/3 from ag- 
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3 Forward- Backward Asymmetry 



The formalism used in the last section to investigate the total cross section for tt produc- 
tion addressed a very suitable observable to determine the mass of the top quark. But 
since the total cross section at threshold is also significantly influenced by other quan- 
tities like the strong coupling constant and the top decay-width, it is necessary to have 
some independent observables. It was first proposed in that the forward backward 
(FB) asymmetry would provide such a quantity. It has the additional advantage that it is 
independent of the absolute normalization of the cross section measurement. Thus it can 
be determined experimentally with high accuracy. However, treating this problem in a 
nonrelativistic context, one encounters unphysical divergences which made necessary the 
introduction of an unnatural cut-off. In this section we will show a systematic way how 
to avoid such divergencies. 

The FB asymmetry is defined as the relative difference of the cross section for particles 
produced in the forward and in the backward direction, with respect to the direction of 
the e~ beam. 



ApB ■= = (36) 

<7 tot 0" tot 

Since the top quarks will be identified by their main decay products W and b let us 
consider the cross section for the process e'^e^ —>■ ti ^ hW^hW~ . 

da^+^-^bw+bw- = ^\Ma + Mv\'^d^w+bd^w--b{'^^YS{ti + t2 - h - k2 - h - 62) (37) 

with the phase space factor 

d<^w^b = (0i(2^)0(^io)<5(^^ - M'){2n)e{h,)6{bl) (38) 

and a similar one for the decay of the antitop. Production of W^hW~h via other channels 
can be treated as in 11121 . 



8 




The FB asymmetry originates in the interference terms of the vector couphng (My) 
and the axial couphng (M4). This is depicted in figure (|l|) where V indicates a vector and 
A an axial vector coupling, the star denotes complex conjugation. The matrix element 
M is essentially the vertex '-yti and Zti. Near threshold the perturbative treatment of 
the Coulomb interaction is no longer valid. Instead one has to include the whole rung of 
Coulomb interactions, and even worse, it is necessary to include also the leading logarith- 
mic contributions from the gluonic self energy corrections. We can sum the relevant set 
of graphs by means of the equation 

= - iKDT^ (39) 

for the vertex function T^. D contains the unconnected fermion propagators and K rep- 
resents the sum of all 2pi graphs. Comparison with the BS-Equation for the Feynman 
amplitudes 

iC = -D + DKG (40) 
shows that eq. ( |5UD has the solution 

= -iD-^GY (41) 
= -^D-'G^,Y, (42) 

where it is understood, that the two legs of G on the right hand side are connected with 
the 7-matrices and a momentum integration is performed. 
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We can split the matrix element into 

M\ = mim2DT\ (43) 
M'y = mim2Dry 

where 

mi = e'^;){k,)u,{h){-t " 1^337"^-) (44) 

m, = e*y\h){-i ^ V:,,YP-)vAb,). (45) 
V ^ sm Uw 

Near threshold the propagators in D and the Green function G should be dominated by 
small momenta and therefore it should be possible to replace them by their nonrelativistic 
approximations. However, for D the approximation to 0{a) is needed for the axial vertex, 
since the zero order propagators would give zero due to A~757*A+ = 0. 

D ^ (46) 
This gives rise to the effective nonrelativistic axial vertex 

k 

75y-^[y,7l75A-. (48) 
We note already here that the term ^ introduces an additional power in p which is not 

k 

present in the relativistic propagator. Instead the factor ^ in (^) would be replaced by 
^ which is finite for p ^ oo. Thus a divergence of (|^) at |p| — oo is really an artifact 
of the nonrelativistic approximation. 

Using the fact that the heavy quarks are on shell up to O(a^) we can write 



j d^w+bY.\^^^ = -2ImE(m2)A+ + 0{a^) (49) 

(T,A 



in agreement with [p!3| . 

From the modulus squared of the propagators we obtain a factor 

dpo . 1 ,2 ^ 2 
27r'p2_^2l - r|2cu|2' 



(50) 



uj = (p'^ — mE — imT) (51) 
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Collecting everything from above, performing the trace in a frame where the leptons 
move along the z-axis we get for the FB-asymmetry 



ISr f Sq q^^^ 

vr^r^^ JO J {It^}^ m 



where 



= a(^)6(^) + a(^)6(^) (53) 
d^y = a(^)a(^) + (54) 
J d^FB = jj d^{j^ de - j"^ de)sm9 (55) 



The SM values for a(^),6(^) are given in (|^ and (H). We introduced an intermediate 



cut-off A in (|5^ ) since the p integration is logarithmically divergent. It will be possible to 



remove it in the following. The Green function G{p) in (|52D is defined by 

/d^q Att f°° 

-— rG*(p, g) = — / drr sin prGi=o{r,0). (56) 
[Zn)-^ p Jo 

Gi=o{r,r') denotes the S-wave Green function in configuration space, eq. (pQ]). Let us now 
investigate the term 



/^G(p, q)^— = -f d'xe-^^^^G{x, y)\.^, (57) 

Using the representation (|1^) one can show that only P- waves (/ = 1, m = 0) contribute 
to the sum. Thus we need to evaluate eq. ( P!U[ ) for Z = 1. The regular and singular solutions 
behave at the origin as 



r 

respectively. We now construct a singular solution with c> = 1 out of two arbitrary 
solutions Upi and Up2 ■ 

g^=\r) = ap[up2{r) - BpUpiir)] (59) 
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by means of 



a. 



r->0 

Bp is determined by tlie requirement of regularity at infinity 



lirnr ^[up2{r) - BpUpi{r)] \ (60) 



= l,m (61) 



Then the condition ( pS]) demands c< = m/3. Performing the different at ion in the z- 



direction in the limit ?/ ^ in (|57|) leads to 



(i^g ^q^^^ COS 9 /■°° , , ,=1, , 

TG[p,q) = — — / ar[ pcospr)gy [r). [bZ) 



(27r)'^ ' m Jo r 

While the expressions (|56|) and are well defined, it has been observed in B and 



IJ] that (|52D is logarithmically divergent. Thus cutoffs (A) have been introduced to 
make numerical predictions possible. Since in our present case the divergence is only 
logarithmical, this leads to phenomenologically reasonable results. But the explanations 
given in these references are different and so are the cut-offs. It may, however, happen that 
for the comparison with experimental data with some cuts a reintroduction of some kind 
of cut-off will be necessary |^, but we feel saver by first giving a satisfactory theoretically 
prediction and leaving this possibility open to the experimentalists. 

In any case, the introduction of these cut-offs leads to some numerical uncertainty in 
the result, and is clearly unsatisfactory from the theoretical point of view. Furthermore 
we will see in the next section, that a similar UV-divergence occurs in the context of the 
stop-anti-stop production. Since this divergence is linear, a better understanding of the 
origin of this kind of divergencies is necessary to give a quantitative prediction. 

The key point of our solution of the problem is to go back to the perturbative sum of 
the graphs for apB- This is illustrated in figure (p. 

Consider now the leading (tree) contribution in the nonrelativistic approximation (e.g. 
the first graph in fig. |I|). In this case both S- and P-wave Green functions are replaced by 

{2nf6{p-p') 

G{p,P )= ^ • (63) 

m 



^The author is grateful to M.Jezabek for this information 
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This leads to the logarithmically divergent expression (for A ^ oo) 

(i),nr- 9r f cosO 

I m I 

Since the remaining graphs in fig. |l| by power counting can be shown to give finite results 
we conclude that the nonrelativistic approximation was not valid in the tree graph due to 
the extra power in p from the axial vertex. Furthermore since (|6^) is the only divergent 
contribution in the (infinite, but assumed to be convergent) sum in fig. ^ representing 
G we conclude that the divergent part of (|5^) is entirely contained in the first (free) 
contribution (|6^. We could now remove this divergence by a replacement p/m ^ p/Ep, as 
indicated above. This, however, does not respect another qualitative difference between the 
exact tree contribution and its nonrelativistic approximation. Namely in the relativistic 
calculation the phase space is cut off when the momentum squared of one quark falls 
below the invariant mass of the decay products |§. 

Therefore, we return, instead, to the relativistic expression for the tree contribution. 
Since we are dealing with a decaying particle, we will have to include the self energy con- 
tribution due to the decay in the propagator. This will in general lead to gauge dependent 
results, but we can use the constant on-shell width to obtain the leading contribution in 
the weak coupling . 

This makes the quantities we are considering finite (for F ^ 0) and gauge independent. 
The remaining higher order contributions (e.g. gauge dependent if one considers only one 
diagram) can be calculated perturbatively. Therefore, to leading order it seems reasonable 
to take the constant width approximation for the top propagator. 

Since we are only investigating processes with ti intermediate states which will give 
the main contribution to the cross sections, we focus on the graph shown in fig. 0. A 
further advantage of this method is that now the other - non resonant - graphs for the 



process e~^e~ W^bW~b to leading order in the weak coupling [|l^] need only to be 
added to yield the background contribution. 

Performing the straightforward calculation for the relativistic tree contribution with 
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a constant, but non-zero width F, we arrive at 



a 



(1) 



'FB - -^{cAV + CavP) J dfll J rf/i^A(/i?)A(/i^; — J — 

where 



(65) 



^(p^) = ~ir~2 — ^ (66) 

This could have been also obtained by replacing the fermion propagator by 

J 7r[(m2 — /i^)^ + m^r^] — /x^ + ie 
and cutting it to effectively replace the phase space element for a stable quark 



(instable = 7^0(Po)5(p^ - m'^] 



(2vr) 
by 

C^^nnsiafeZe = j^Qipo) A{p^) ■ (68) 

A transparent summary of our approach may be formulated like this: To obtain finite, 
gauge independent results to the desired order we simply propose to replace the nonrela- 
tivistic tree contribution by the relativistic one and leave the (finite) rest unchanged: 

(^FB = Cr?B + lin^ {(^FB,A - CTfbT) (69) 

A— ♦oo ' 

In fig. ^the numerical results are shown for the different contributions to 0"^?^ for ^ = 300 
GeV. This can be estimated to give a result for lim\^oo{cFB,A — o'fs'a) ^J^^S oiilj 2% 
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Fig. 3: Different Contributions to ufb 
dashed: (Jfb,k, dottedicx^^Y' dashed-dotted: 



E/GeV 



below the final answer. The results of the purely nonrelativistic calculations are compared 
to our approach in fig. H for nit = 180 and an electron polarization of 0.6. We also included 



the hard corrections as given in [14, 16 1. 

We believe that our approach has the decisive advantage that it can be based upon to 
the original set of graphs to be considered and it is clear which graph has been calculated 
to which accuracy. Thus at least in principle a systematic improvement is possible. The 
present approach should also be applicable to the 0{as) final state corrections calculated 
in PI. 
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4 Axial Contribution to the tt Total Cross Section 
and the Production of Stop-Antistop near Thresh- 
old 

The total cross section for the production of P-wave states near threshold gives the leading 
term for the production of stop-antistop {tt) and a contribution of to the total cross 

section for ti. Only qualitative investigations on the basis of the Coulomb Green function 
have been published up to now [|1^, [1^. We will apply the method of replacing divergent 



nonrelativistic graphs, implicit in the Green function approach, by the relativistic, finite 
ones to get a reliable quantitative result also for this case. 

Considering the Coulomb Green function as a first approximation it was observed 
in |jl8[ that the total cross section for a stop-antistop pair near threshold develops an 
unphysical linear divergence due to the nonrelativistic approximation. Even a relativistic 
calculation of the imaginary part of the one loop contribution using the commonly used 
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propagator — + imT) remains logarithmically divergent. This is due to the fact 

that this propagator does not fulfill the requirements of (local) quantum field theory as 
does the expression 

I , (70) 

which is in agreement with the Lehmann representation of the full propagator. The tree 
contribution (denoted cxi) using (^) reads 



whereas we have for the nonrelativistic tree contribution 

3r 

However, a single subtraction of the tree graph is not sufficient. One should also subtract 
the logarithmically divergent one-gluon exchange term. But since we replaced the gluon 
propagator by the resummed one we should also calculate the relativistic graph with a 
resummed propagator to get the correct behavior at infinity. Unfortunately this leads to 
renormalon ambiguities. Therefore, we conclude that the uncertainties in choosing the 
right cut-off are of the same order of magnitude as higher QCD corrections, and we 
circumvent these difficulties for the time being by keeping the cut-off in the logarithmic 
divergent terms, choosing its value A = Aq = m from the observation that this cut-off 
would have given the correct answer in the case of the FB-asymmetry discussed in the 
last section. 

The same is true for the pure axial contribution to the total cross section for ti production 
since each axial vertex contributes in the nonrelativistic limit an extra power in p as 
explained in the last section. The only change is to replace the factor 

fs--=\{cs + dsp) (73) 

by 

fr.= Q{cA + dAp). (74) 
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Fig. 5: Different Contributions to a^-^ 
dashed: (T^°, dotted: erf dashed-dotted: ai^rei 



Using the methods of the preceding sections we obtain within the Green function approach 



cr 



(75) 



(27r)3 m 

where fx denotes either fs for the scalar or // for the axial fermionic cross section. The 
different contributions to the cross section for the scalar case 



(76) 



are depicted in fig. || for cos^ 9t = 0.5, m = 180, T = Ttop- The net result is shown in fig. |] 
for the axial contribution to the ti cross section. 
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Fig. 6: Drell-ratio (R= caIo-^p) for the production of tt 



The full 0{a) corrections are taken into account by including the hard corrections 
by a factor (1 — pO[ and a 10% reduction of the decay width ||21| in addition to the 
aforementioned potential. We conclude that the axial contribution gives a sizeable effect 
(of a few percent ) to atot beginning at ~ 5GeV above threshold. However, it is small 
enough in order to hide the uncertainty in A with respect to 0{af) corrections to the 
vector contribution to atot- 

In fig. we compare the cross sections for the production of tt near threshold for two 
different values of the decay constant ( cos"^ Of = 0.5, m = 180). 
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Fig. 7: dashed curve: F = l.SGeV, full line: F = O.SGeV 
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It is interesting to note that in contrast to the S-wave case the 0{as) corrections 
to the potential give a relatively small effect. This raises the hope that in the P-wave 
case fixed order perturbation theory works better as it does for S-waves and thus our 
subtraction procedure is stricly applicable also for the one gluon exchange term without 
getting troubles from the resummed gluon propagator. In fig. |] and ^ we compare our 
result for the tt cross section with that for a pure Coulomb potential, both computed with 
our subtraction method. 

R 





0.08 












0.06 












0.04 












0.02 , 










-4 


-2 


2 


4 


6 


3 10 



E/GeV 



Fig. 8: dashed curve: Coulomb potential; full line: QCD-potential; both with 

r = l.SGeV 
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Fig. 9: dashed curve: Coulomb potential; full line QCD-potential; both with F = O.SGeV 
A further difference is that the lowest lying resonance peak (the 2P-level) is enhanced 
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by the QCD potential and not suppressed as is the IS-peak. This is in agreement with 
the calculation of the wave function corrections [p2 |. 



5 Conclusion 

We have calculated the forward-backward asymmetry and the pure axial contribution 
to the total cross section for the tt system. Furthermore the leading contribution to the 
total cross section for heavy scalars is determined in this region. Since the necessity 
of the inclusion of an infinte sum of Coulomb interactions in practise enforces the use 
of a nonrelativistic Schrodinger equation, unphysical divergencies are produced in the 
nonrelativistic limit. These are removed by showing that they originate in the leading 
nonrelativistic contributions which may be subtracted and subsequently replaced by the 
relativistic expressions which lead to finite results. Thus we are able to combine both 
the advantages of a nonrelativistic calculation near theshold with the one of a strictly 
perturbative treatment. 

Acknowledgement: This work has been supported by the Austrian Science Founda- 
tion (FWF), project P10063-PHY within the framework of the EEC- Program "Human 
Capital and Mobility", Network "Physics at High Energy Colliders", contract CHRX- 
CT93-0357 (DG 12 COMA). I am grateful to Y.Sumino for sending me the essential parts 
of his Fortran code for the numerical evaluation of the Green function. I am also grateful 
to M.Jezabek, V.A. Khoze and J.H. Kiihn for instructive discussions and to W. Kummer 
for a careful reading of the manuscript and many useful suggestions. 



21 



References 

[1] F.Abe et.al. (CDF coll.):" Observation of Top Quark Production in pp Collisions", 
FERMILAB-PUB-95/022-E. 

[2] S.Abachi et.al. (DO coll.):" Observation of the Top Quark" , FERMILAB-PUB-95/354- 
E, liep-ex-9503003. 

[3] G: Altarelli:" Phenomenology in and beyond the Standard Model", CERN-TH/95-3, 
and references therein. 

[4] V.S. Fadin, V.A. Khoze, Sov.J. Nucl.Phys.48 (1988) 309. 

[5] Y. Sumino, K. Fujii, H. Hagiwara, H. Murayama, C.-K. Ng, 
Phys.Rev. D47 (1993) 56. 

[6] M. Jezabek, J.H. Kiihn, T.Teubner, Z.Phys.C56 (1992) 653; 
M. Jezabek, T.Teubner, Z.Phys.C59 (1993) 669. 

[7] M.J. Strassler, M.E. Peskin, Phys.Rev.D43 (1991) 1500. 

[8] "Gluon Radiation Off Scalar Stop Particles", W.Beenaker, R.Hopker, P.M.Zerwas, 
DESY 94-235 

[9] H. Murayama, Y. Sumino, Phys.Rev.D 47 (1993) 82 

[10] A.Billoire, Phys.Lett.B 92 (1980) 343. 

[11] W.Kummer, W.Mddritsch, Z.Phys. C66 (1995) 225. 

[12] A. Ballestrero, E. Maina, S. Moretti, Phys.Lett. B 333 (1994) 434 

[13] W. Modritsch, W.Kummer, Nucl. Phys. B430 (1994) 3; W.Kummer, W.Modritsch, 
Phys.Lett B349 (1995) 525. 

[14] R.Harlander,M. Jezabek, J.H. Kiihn, T.Teubner, "Polarization in Top Quark Pair Pro- 
duction near Threshold", |hep-ph/9411395| , TTP94-28. 



[15] M. Veltman, Physica 29 (1963) 186, H. Veltman, Z.Phys.C 62 (1994) 35. 
[16] J.H. Kiihn, P.W. Zerwas, Phys.Rep.167 (1988) 321. 
[17] Y.Sumino, Ph.D.Thesis Univ. of Tokyo, 1993, UT-655 

22 



[18] I.I.Bigi,V.S.Fadin,V.Khoze,Nucl. Phys. B377 (1992) 461. 

[19] V.S Fadin, V.A. Khoze, Sov.J. Nucl.Phys.53 (1991) 692; (Yad.Fiz.53 (1991) 1118). 

[20] S.Giisken, J.H.Kiihn, P.M.Zerwas, Phys.Lett. B155 (1985) 185; 
J.Jersak, E.Laermann, P.M.Zerwas, Phys.Rev.D25 (1982) 1218. 

[21] M.Jezabek, J.H.Kiihn, NucLPhys. B314 (1989) 1. 

[22] W.Kummer, G.Wirthumer, NucLPhys. B185 (1981) 41, E:B194(1982) 546 



23 



